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We present 17 experimentally confirmed Feshbach resonances in optically trapped 85 Rb. Seven 
of the resonances are in the ground-state channel (/, rrif) = (2, +2) + (2,-1-2), and nine are in the 
excited-state channel (2, —2) + (2, —2). We find a wide resonance at high field in each of the two 
channels, offering new possibilities for the formation of larger 85 Rb condensates and studies of few- 
body physics. A detailed coupled-channels analysis is presented to characterize the resonances, 
and also provides an understanding of the inelastic losses observed in the excited-state channel. In 
addition we have confirmed the existence of one narrow resonance in a (2, +2) + (3, +3) spin mixture. 



I. INTRODUCTION 

The creation of ultracold molecules is currently of great 
interest. They offer a wide range of applications includ- 
ing: studies of few-body quantum physics, high preci- 
sion spectroscopy, quantum simulators for many-body 
phenomena and controlled chemistry [TJ [2]. Ultracold 
molecules have a far richer substructure than atoms and 
so molecular condensates with tunable interactions offer 
unique levels of control over collision properties [3]. One 
route to ultracold molecules is through the association 
of two ultracold atoms into a weakly bound molecule 
[3]. The energy of a bound molecular state is tuned 
adiabatically through an avoided crossing with the en- 
ergy of the separated atomic states [3J , forming a weakly 
bound molecule. The molecules can then be transferred 
into their ro-vibrational ground state by stimulated Ra- 
man adiabatic passage (STIRAP). This method has been 
used effectively in several systems to create ultracold 
molecules [5HZ]. 

85 Rb is a promising species for ultracold atomic gas 
experiments, though it has often been overlooked due 
to the challenges of forming a Bose-Einstein condensate 
(BEC) [HI IS]- Our recent work shows the benefits of 
85 Rb for RbCs production [ID]. However, for these ex- 
periments a full understanding of the scattering behav- 
ior of 85 Rb is required. Most previous work on 85 Rb 
has focused on the wide resonance near 155 G in the 
(f,m f ) = (2, -2) + (2, -2) channel [IT]. This resonance 
is suitable for experiments that require precise tuning of 
the scattering length and has been used extensively in 
studies of condensate collapse [DJ [T2"Hh"i] , the formation 
of bright matter wave solitons [T5] , and few-body physics 
|16| . Further work using 85 Rb includes spectroscopic 
studies of photo-association [17] [18] and measurements 
of inelastic collision rates [TD1 |2"D] , molecular binding en- 
ergies [UJ , molecule formation [2"2"H2"rJ] and Efimov states 
[T4l [TBI [27] . Despite extensive work in this region of the 
(2, —2) + (2, —2) channel, there appears to have been lit- 
tle theoretical or experimental work on the ground state, 



or at higher field. 

In this paper we reveal the rich Feshbach structure of 
85 Rb. We use coupled-channels calculations to predict 
Feshbach resonances in both the (2, —2) + (2, —2) chan- 
nel (designated ee), and (2, +2) + (2, +2) channel (desig- 
nated aa) and confirm 16 of them experimentally. In ad- 
dition we identify a resonance in the mixed spin channel 
(2, +2) + (3, +3). The structure of the paper is as follow: 
Section II describes the theory and calculations; Section 
III describes the experimental setup and methodology; 
Section IV describes the results, including an outlook on 
future research prospects. 

II. THEORY 

The collision Hamiltonian for a pair of alkali-metal 
atoms is 




+ H 1 +H 2 + V(r), (1) 



where r is the internuclear distance, fi is the reduced 
mass, L is the rotational angular momentum operator 
and V is the interaction operator. Hi and H 2 are the 
monomer Hamiltonians of the free atoms, 

Hi = (iH ■ §i + (g c fi B Si Z + g n ^BHz)B, (2) 

where Q is the hyperfine coupling constant of atom i, g c 
and g n are the electron and nuclear g-factors, s and i 
are the electron and nuclear spin operators and B is the 
magnetic field. 

The calculations in the present paper are carried out 
in two different basis sets: a fully decoupled basis set 

IsRbmsRb) \iRbmiBb) |sRbW sRb ) |iRb"T- iRb ) \LM L ) , 

and a partly coupled basis set 

\fa,m fia ) \ f b ,m f . b ) \F,Mf) \L,m L ) . 
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FIG. 1. (Color online) Upper panel: The calculated s-wave scattering length in the aa channel of 85 Rb2, with resonances 
marked with lines whose color depends on their Mp value (see legend); the length of each line is proportional to the logarithm 
of the width of the resonance. Lower panel: The energies of weakly bound molecular states, relative to the aa threshold, 
(2, +2) + (2, +2) channel. All calculations in this figure are for Mtot = 4, corresponding to s-wave scattering in the aa channel. 



The two basis sets give identical bound-state energies and 
scattering properties, but different views of the bound- 
state wavef unctions. In both cases the basis sets are 
symmetrized to take account of identical particle sym- 
metry. The resulting coupled equations are diagonal in 
the total projection number M tot = Mp + Ml, where 
M F = m/ l0 + m/,& = m sRba + TOiRb a + m sRbb + m iRbb . 
The basis sets used include all functions with L = and 2 
for the required value of M tot , which for s-wave scattering 
is equal to my -a + m f,b for the incoming channel. 

The coupled-channel scattering calculations are per- 
formed using the MOLSCAT program [55], as modified 
to handle collisions in an external field [35]. Calcu- 
lations are carried out with a fixed-step log-derivative 
propagator |30j from 0.3 nm to 2.1 nm and a variable- 
step Airy propagator [31 from 2.1 nm to 1,500 nm. 
The wavefunctions are matched to their long range so- 
lutions, the Ricatti-Bessel functions, at 1,500 nm to 
find the S-matrix elements. The s-wave (L = 0) scat- 
tering length a(k) is then obtained from the identity 
a(k) = (ik)- 1 ^ - 5 00 )/(l + Sao) [32], where S 00 is the 
diagonal S-matrix element in the incoming channel and 
k is the wavevector. The bound-state calculations use 
the BOUND [33J and FIELD [3J] packages, which locate 
bound states using as described in Ref [35] • BOUND 



and FIELD use propagator methods, without radial basis 
sets. The calculations allow the assignment of quantum 
numbers to the states responsible for resonances in the 
scattering length. 

The scattering and bound-state calculations are car- 
ried out using the potential curves and magnetic dipole 
coupling function from Ref. [36] . The potentials were 
obtained by fitting to spectroscopic data on both the 
singlet [37] and triplet states of 87 Rb2 and the triplet 
state of 85 Rb2, together with several Feshbach resonances 
in 87 Rb 2 , 87 Rb 85 Rb and 85 Rb 2 . The singlet and triplet 
scattering lengths for 85 Rb on the potentials of ref. [3S] 
are as — 2735 clq and ap = —386 do respectively. 

The calculated s-wave scattering length for the aa 
channel is shown in the top panel of Figure [I] and the 
binding energies of the near-threshold molecular states 
responsible for the resonances are shown in the lower 
panel. The resonance positions are given in Table [I] 
along with their widths A as defined by local fits to 
the standard formula a(B) = a bg [1 — A/(S — Bo)] |38j . 
where a bg is the background scattering length, A is the 
width, and Bq is the position of the pole in the scattering 
length. Quantum numbers were assigned by carrying out 
approximate calculations with either Mp or F and Mp 
restricted to specific values. For a homonuclear diatomic 
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FIG. 2. (Color online) Upper panels: The real and imaginary parts of the s-wave scattering length in the ee channel of 85 Rb2. 
Each resonance is indicated by a colored vertical line that indicates its Mf value (see legend); for pole-like resonances, the 
length of the line is proportional to the logarithm of the width of the resonance. Inelastically dominated resonances are not 
always evident in a but appear as peaks in /3 and are indicated by dashed vertical lines. Lower panel: The energies of weakly 
bound molecular states, relative to the ee threshold, (2, —2) + (2, —2) channel. Only states with no continuum interference 
(Mf < —4) are shown in the bound-state map, but all resonances are included in the scattering length. All calculations in this 
figure are for M tot — —4, corresponding s-wave scattering in ee channel. 



molecule, F is a nearly good quantum number in the low- 
field region where the free-atom energies vary linearly 
with B. Figure [l] shows one wide resonance near 851 G 
(A = 1.2 G) that offers attractive possibilities for pre- 
cise tuning of the scattering length, and many narrower 
resonances that may be useful for molecule formation. 

For an excited-state channel, where inelastic scatter- 
ing can occur, the scattering length a(B) is complex, 
a(B) = a(B) — i/3(B). The two-body inelastic loss rate 
is proportional to /3(B). The upper panels of Figure [2] 
show the real and imaginary parts of a(B) for s-wave 
collisions in the ee channel. In this case the inelastic col- 
lisions produce atoms in lower magnetic sublevels, with 
m/ jQ and/or m/,& > —2. The lower panel shows the cor- 
responding molecular bound states for Mp = —4, —5 and 
—6, obtained from calculations with Mp fixed. We also 
carried out calculations of the quasibound states with 
Mp = — 2 and —3 near the ee threshold in order to iden- 
tify the states responsible for the remaining resonances. 
These calculations use the FIELD program with propa- 
gation to reduced distances around 100 nm in order to 
reduce interference from continuum states. 

In the presence of inelastic scattering, a(B) does not 
show actual poles at resonance [32]. If the background 



inelastic scattering is negligible, the real part a(B) shows 
an oscillation of amplitude a rcs , while the imaginary 
part shows a peak of height a rcs . The resonant scat- 
tering length a ros is determined by the ratio of the cou- 
plings from the quasibound state responsible for the res- 
onance to the incoming and inelastic channels [32]. If 
there is significant background scattering, then there is a 
more complicated asymmetric lineshape that may show 
a substantial dip in the inelastic scattering near reso- 
nance [3JJ] ■ Figure [2] shows resonances of all these dif- 
ferent types: the resonances due to bound states with 
Mp = mf,a + m f.b = — 4, —5 and —6 are pole-like, with 
values of at least a res > 20 a and with most a les > 1000 
aQ. These resonances produce pronounced features in 
a(B) and sharp peaks in /3(B), off scale in Figure [2] By 
contrast, resonances due to states with Mp = —2 and 
—3 show much weaker features with a rcs < 15 a and 
some lower than 0.01 arj. These are barely perceptible 
in a(B) on the scale of Figure [5] and produce broader, 
weaker peaks in P(B). The distinction occurs because 
all the inelastic channels have Mp > —4: bound states 
with Mp = — 2 and —3 are generally more strongly cou- 
pled to inelastic channels with the same Mp than to the 
incoming channel with Mp = —4, whereas the reverse is 



4 



true for bound states with Mp = —4, —5 and —6. Many 
of the features show quite pronounced asymmetry in the 
shape of the inelastic peaks. All of the resonances with 
a ros > 1.0 ao are listed in Table [TT| along with their widths 
and approximate a rcs values. 

We have also investigated the scattering length for 
mixed spin channels with a view to identify broad res- 
onances suitable for manipulating interactions. Most 
channels exhibit strong inelastic decay with measured 
trap lifetimes of ~ 100 ms. However, the (2, +2) + (3, +3) 
channel is immune to inelastic spin exchange collisions, 
resulting in trap lifetimes of ~ 5 s. The scattering length 
in the mixed spin channel, (2, +2) + (3, +3) shows two 
pole-like resonances at 818.8 G and 909.9 G, both with 
widths of 2 mG, and a rcs — 1600 ao and 800 ao respec- 
tively. 




Magnetic Field (G) Magnetic Field (G) 

FIG. 3. Narrow resonances in 85 Rb, with fitted width 5 < 
0.2 G, observed as loss features in the atom number, (a) A 
resonance in the (2, —2) state at 248.64(1) G. (b) A resonance 
in the (2, +2) state at 770.81(1) G. The error bars show the 
standard deviation for multiple control shots at specific mag- 
netic fields. 



III. EXPERIMENT 

The details of our apparatus and cooling scheme are 
presented in Refs. HOI SI] and are only briefly re- 
counted here. Ultracold samples of 85 Rb are collected in a 
magneto-optical trap before being optically pumped into 
the (2, —2) state and loaded into a magnetic quadrupole 
trap. Forced RF evaporation cools the 85 Rb atoms to 
50 /iK where further efficient evaporation is impeded by 
Majorana losses [32]. The atoms are then transferred 
into a crossed dipole trap derived from a single-mode 
1550 nm, 30W fibre laser. When loading, the power in 
each beam is set to 4 W, creating a trap 100 /xK deep 
with radial and axial trap frequencies of 455 Hz and 
90 Hz respectively. After loading, when performing Fes- 
hbach spectroscopy in the absolute internal ground state, 
the 85 Rb atoms are transferred into the |2, +2) state by 
RF adiabatic passage [13] ■ A vertical magnetic field gra- 
dient of 21.2 G/cm is then applied, slightly below the 
22.4 G/cm necessary to levitate 85 Rb. In contrast, when 
working with the (2,-2) state, no magnetic field gradient 
is applied and the atoms are confined in a purely optical 
potential. 

A typical experiment begins with 6.0(3) x 10 5 85 Rb 
atoms at 10.2(1) /iK confined in the dipole trap in either 
(2, +2) or (2,-2). To perform Feshbach spectroscopy, 
the magnetic field is switched to a specific value in the 
range to 1000 G. Evaporative cooling is then performed 
by reducing the dipole beam powers by a factor of 4 over 
2 s. The atomic sample is then held for 1 s in this fi- 
nal potential. Resonant absorption imaging is used to 
probe the atoms after each experimental cycle. Feshbach 
resonances are identified by examining the variation in 
the atom number and temperature as a function of the 
magnetic field. The magnetic field is calibrated using 
microwave spectroscopy between the hyperfine states of 
85 Rb. These measurements reveal our field stability to 
be 0.1 G for the range to 400 G and 0.5 G for the range 
400 to 1000 G. 

To perform Feshbach spectroscopy on a spin mixture 



of (2, +2) + (3, +3), a cloud of (2, +2) atoms is first cooled 
using the same evaporation sequence as above at 22.6 G. 
To populate the (3, +3) state a microwave pulse is applied 
for 250 ms at 3.0887 GHz producing a mixture containing 
7(1) x 10 4 atoms in each of the (2, +2) and (3, +3) spin 
states. The magnetic field is then switched to a value 
in the range to 1000 G and held for 750 ms. Finally, 
Stern-Gerlach spectroscopy and absorption imaging are 
used to probe both spin states simultaneously. 

IV. EXPERIMENTAL RESULTS 

We have observed 7 resonances in the aa channel and 
9 resonances in the ee channel. The observed and pre- 
dicted resonance positions and widths are listed in Tables 
[I] and [H] and show good agreement between experiment 
and theory. In the ground state, all the widest calculated 
resonances are seen experimentally, with the exception 
of the two high-field resonances where the experimental 
field is less stable. In the excited state, all but two of the 
predicted pole-like resonances are seen, together with two 
of the inelastically dominated features. 

Figure [3] shows fine scans of the atom number for two 
of the narrow resonances, one in each of the aa and ee 
channels. Such resonances produce sharp drops in atom 
number. The experimental positions and widths of these 
resonances are determined by fitting a Lorentzian, with 
width 5, to the data points. It should be noted that the 
experimental and theoretical widths are entirely differ- 
ent quantities for narrow resonances, and should not be 
compared. 

There are three resonances with widths greater than 
1 G. Figure [4] shows a fine scan across the resonances 
near 530 G in the ee channel, and 850 G in the aa chan- 
nel. In these cases the atom number shows both a peak 
and a trough. The trough (loss maximum) again corre- 
sponds to the resonance position, while the peak (loss 
minimum) occurs near the zero-crossing of the scattering 
length. The three wide resonances are several orders of 
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FIG. 4. Broad resonances in 85 Rb, with width A > 1 G, 
observed as features in the atom number, (a) A resonance in 
the (2, -2) state at 532.3(3) G. (b) A resonance in the (2, +2) 
state at 852.3(3) G. The experimental widths are determined 
by the difference between the positions of the minima and 
maxima in the atom number marked with solid lines in both 
plots. The error bars show the standard deviation for multiple 
control shots at specific magnetic fields. 
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TABLE I. Location and assignment of Feshbach resonances 
for 85 Rb2 in the aa channel in the field range between and 
1000 G. All quantum numbers in the table refer to the molec- 
ular states. The experimental errors shown are statistical 
uncertainties resulting from the fits as described in the text. 
The experimental width labeled ' was determined from the 
difference between the minima and maxima in the measured 
atom number. Additional systematic uncertainties of 0.1 G 
and 0.5 G apply to resonance positions in the field ranges 
to 400 G and 400 to 1000 G respectively. 



magnitude wider than any of the other resonances seen 
and provide valuable control over the scattering length. 
Note that our measurement of the position of the well- 
known 155 G resonance in the ce channel is not as accu- 
rate as the determination from bound-state spectroscopy 
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FIG. 5. The two inelastically dominated features observed in 
the (2, —2) state. Figures (a) and (b) show the atom num- 
ber while (c) and (d) show temperature. Error bars show the 
standard deviation for multiple control shots at specific mag- 
netic fields, (e) and (f) show the calculated rate coefficient 
for two-body loss K\ oss . 



experimentally are those with the largest a rcs values. The 
number of atoms in the trap decreases around these res- 
onances due to an increase in the 2-body loss rate, as 
shown in Figure ^ a) and (b). The inelastic collisions 
also lead to an increase in temperature, as shown in Fig- 
ure [5] (c) and (d) . The rate coefficient for 2-body losses 
due to inelastic collisions from a channel n is 



2h 



-9nP{B), 



(3) 



The two inelastically dominated features that are seen 



where /3(B) is the imaginary part of the scattering length 
and g n = 2 for a thermal cloud of identical bosons [3]. 
The calculated rate coefficients for the two resonances are 
shown in panels (e) and (f); they peak around 1 x 1CP 11 
cm 3 /s, which is an order of magnitude higher than for 
any of the other inelastically dominated features. 

We have also measured one resonance in the (2, +2) + 
(3, +3) mixed spin channel. The experimental results are 
presented in Figure [6] where a loss feature in the 85 Rb 
(2, +2) number reveals the location of the resonance. A 
Lorentzian fit gives a resonance position of 817.45(5) G 
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FIG. 7. The rate coefficient for 2-body loss K\ OBS (red solid 
lines), which is proportional to the imaginary part of the scat- 
tering length, and the corresponding real part of the scatter- 
ing length (dashed lines) for the two resonances with A > 1 G 
in the ee channel. Note the dip in -K"i oss on the high-field side 
the 532 G resonance. 



TABLE II. Location and assignment of Feshbach resonances 
for 85 Rb2 in the ee channel in the field range between and 
1000 G. All resonances shown satisfy a ICS > 1 ao- AH quan- 
tum numbers in the table refer to the molecular states. The 
experimental errors shown are statistical uncertainties result- 
ing from the fits as described in the text. The experimental 
width labeled ' was determined from the difference between 
the minima and maxima in the measured atom number. Ad- 
ditional systematic uncertainties of 0.1 G and 0.5 G apply to 
resonance positions in the field ranges to 400 G and 400 to 
1000 G respectively. The resonances near 155 G and 220 G 
have been measured previously |21l 144] . 
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FIG. 6. A resonance measured between the (2, +2) and 
(3, +3) spin states in 85 Rb at 817.45(5) G. On resonance the 
increased inelastic collision rate in the mixture results in a 
loss feature in the (2, +2) atom number as a function of mag- 
netic field. The error bars show the standard deviation for 
multiple control shots at a specific magnetic field. 

and an experimental width of 0.031(1) G, which may be 
compared with the predicted position 818.8 G. 



V. CONCLUSION 

A detailed understanding of the two-body scattering 
behavior is essential for understanding many phenom- 



ena in ultracold gases. These include studies of molecule 
formation [22H2S], Efimov states and their universality 
[HI ITrJI |2"71 25], dimer collisions and few-body physics 
[46] , BEC production [9] 07] , controlled condensate col- 
lapse [T2HT4] . and the formation of bright matter wave 
solitons [TS]. The scattering properties of many alkali- 
metal atoms have been documented in the literature [3J. 
However, nearly all previous work on 85 Rb has focussed 
on a single broad resonance near 155 G. This paper re- 
dresses this balance by presenting a detailed study of the 
scattering properties of 85 Rb, revealing additional broad 
resonances and numerous unreported narrow resonances 
in both the ee and aa channels. As has been the case for 
other alkali-metal atoms, this work will facilitate many 
future studies using 85 Rb. 

We have recently explored interspecies Feshbach reso- 
nances in mixtures of 85 Rb and 133 Cs [10], as a key step 
towards the production of 85 Rb 133 Cs molecules. The im- 
proved understanding of the collisional behavior of 85 Rb 
resulting from the present work is essential for the pro- 
duction of the high phase-space density mixtures required 
for efficient molecule formation. In particular, the two 
previously unmeasured broad resonances presented here 
offer new magnetic field regions for evaporative cooling. 
The elastic to inelastic collision ratio in the vicinity of 
these features is potentially more favorable for evapora- 
tive cooling than near the 155 G resonance, where direct 
evaporation of 85 Rb to BEC is possible EE]- Figure 
[7] compares the scattering properties around the 532 G 
resonance with those near the 155 G resonance. The re- 
sults for the 562 G resonance show a pronounced dip in 
the rate coefficient for 2-body loss near 570 G, due to 
interference between the resonant and background con- 
tributions to the inelastic scattering [35] , which offers 
a range of magnetic fields where more efficient cooling 
may be possible. No such dip in the 2-body loss rate 
is present near the 155 G resonance. Alternatively, the 
aa channel offers the prospect of evaporative cooling free 
from two-body loss. Although the background scattering 
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length is moderately large and negative for ground-state 
atoms (see Figure [T]), the broad resonance at 851 G may 
be used to tune the scattering length to modest positive 
values, improving the evaporation efficiency and offer- 
ing the prospect of BEC formation directly in absolute 
ground state. In the future we will investigate evapora- 
tive cooling in these new field regions. We will also use 
the improved knowledge of the scattering of 85 Rb pre- 
sented in this paper, together with similar knowledge for 
133 Cs [49j [50], to devise a route to cooling 85 Rb- 133 Cs 



mixtures to suitable phase-space densities for magneto- 
association using a narrow interspecies resonance [10] . 
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